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Abstract 



We propose a bilocal field theory for mesons in two dimensions obtained as a kind of 
non local bosonization of two dimensional QCD. Its semi-classical expansion is equivalent 
to the 1/N C expansion of QCD. Using an ansatz we reduce the classical equation of motion 
of this theory in the baryon number one sector to a relativistic Hartree equation and solve 
it numerically. This (non topological) soliton is identified with the baryon. 
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Its widely believed that quantum chromo dynamics (QCD) is the theory of strong 
interactions. However the degrees of freedom of its original lagrangian (quarks and gluons) 
do not correspond to the physical particles (hadrons). Besides, many properties of the 
physical particles have an intrinsically non-perturbative character, making difficult to make 
predictions starting from first principles. One possible strategy is to construct effective 
theories valid in the low energy domain. In four dimensions the spontaneous breakdown 
of chiral symmetry leads to the non-linear sigma model that describes the lightest mesons 
(pions), that appear in this context as Goldstone bosons. An interesting fact about this 
model is the presence of soliton solutions (the so called Skyrme term [1] has to be included 
to stabilize these solutions). After the inclusion of the Wess-Zumino term [2], this soliton 
is a fermion, has baryon number one, and consequently, can be identified with the nucleon. 
In this way, although nucleons are not explicit in the lagrangian, their properties can 
be studied, with reasonable agreement with experiment [3]. Arguments using the 1/N C 
expansion [4] also suggest that nucleons are solitons [5]. In some aspects two dimensional 
QCD resembles the four dimensional case (for example, it is confining) but it is much 
simpler, providing a place to test ideas too difficult to apply directly in four dimensions, 
and it has been numerically solved using a discrete light cone quantization [6] . One can take 
an attitude similar to the four dimensional case and consider low energy effective theories 
[7], even though chiral symmetry is not spontaneously broken. But in two dimensions 
one can do better: to write an equivalent lagrangian involving only meson fields without 
relying on low energy approximations. The price to pay for this is that the mesons fields 
are not local. To see this let us start with the QCD action with one flavor 



and choose the gauge A- = 0. Only x is a true degree of freedom. In fact, we can use r\ 
and A + equations of motion 
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to write 77 and A°^_ in terms of \ 

Tfl f 

V(x) = -i— / dyh(x - y)x(y) 

, J (4) 
A a + {x) = 9 - J dyG(x-y):^T a X (y):, 

where 

h(x) = -sgn(x), 

\ (5) 

G{x) = 

The hamiltonian (generator of x + translations) we obtain is 

H = -i— / dx dyh(x - y)xL(x)Xm(y) 

+ — / dxdyG(x - y) : xL( x ) T mnXn(x) : : xl{v)Tp q Xq{v) i 
where m,n . . . are color indices. The normal order is in relation to the vacum defined as 

xCp)|o> = o 

p>0, (7) 

X t (-p)|0> = 0. 



with 



Defining 



*(*) = / ^e~ ipX x(p)- (8) 



M(x+,x ,y ) = — )Xm(^,«/ ) (9) 

we can write the hamiltonian(6), after some rearrangement of the interaction term as 

H = —i ^ - c y dx dyh{x — y)M{x,y) 
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with 



(10) 

dx dyG(x — y)M(x, y)M(y, x), 



2 2 9 2 N C 
\i = m . (11) 
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The bilocal field M(x, y) creates a quark-antiquark pair (a meson) at two causally con- 
nected points so a theory in terms of these objects do not violate causality. The commu- 
tation relation for the fields M(x,y) are those of the infinite dimensional general linear 
group with a central extension 

[M(x, y), M(z, u)} = ±-[6(y - z)(M(x, u) - e(x - u)) - 8(x - u) (M(z, y) - e(z -y))], (12) 

where 

«,) = (13) 

One particular representation of (12) corresponds to two dimensional QCD. In order to 
have a theory with no negative norm states and positive energy this representation should 
be chosen unitary and of heighest weight. We can choose the representation this way 
only when iV c is an integer [8] (from the point of view of meson theory that is the only 
reason this should be so). The hamiltonian (10) and the commutation relations (12) define 
a quantum theory of mesons equivalent to QCD, but written entirely in terms of color 
singlets. It is also Lorentz invariant with the transformation rules 

x + ^Xx + 

x'^X^x- (14) 
M(x,y) -> XM(x,y), 

with A related to the boost velocity v by A = (1 — v)(l — v 2 ) -1 / 2 . Notice that the factor 
1/N C plays the role of Tl. Consequently the classical version of the theory defined above 
corresponds to QCD at N c — > oo. The classical equation of motion for M{x,y) is 

d + M(x, y) = ^- J dzM(x, z)h{z - y) - h(x - z)M(z, y) 

J dzG(z-x)M(x,z)(M(z,y)-e(z-y)) (15) 

-G(z - y)M(z, y)(M(x, z) - e(x - z)) . 
The vacuum expectation value of M(x,y) in QCD is zero. If we expand (15) around 
M(x, y) = and keep up to linear terms we have, in momentum space 

d+ M(p, q )= l -^(---)M(p,q) 

\ q P (16) 
ia 2 N f 

— \ dk G(fc)(sgn(p) M(p -k,q-k)- sgn(q) M(p + k,q + k)). 



Defining 

P-=p-q 

t=jr ( 17 ) 

•qip/PJ) =M(p, q) 

and remembering that the Fourier transform of G(x) is a distribution defined by 



/ 



dpG(p)f(p) = J dp±(f(p)-f(0)). (18) 



we have 
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(19) 



This is just the equation derived by tHooft [9] summing the "rainbow" diagrams whose 
solution gives the meson spectrum (up to order of 1/N C ). As shown by tHooft there are 
no scattering states, implying confinement. 
The baryon number can be writen as 

B = J dxM(x,x). (20) 

Thus the nucleon should be identified with the lowest energy static solution of (15) such 
that B defined by (20) is equal to 1. Due to the complexity of (15) we will use numerical 
analysis and the ansatz 

M(x,y) = ^*(x)^(y) + e(x-y), (21) 



with 



/*W*)|» = 1, (22) 



where \1/ is a classical bosonic variable. With this ansatz the equation of motion (15) 
reduces to an equation for W 

d+V(x) = ^- J dzh(z -x)V(z) +zy J dzG{z - x)\^{z)\H{x) (23) 



Note we have now m instead of the previous renormalized mass \i. A static M(x, y) does 
not imply a static ^(x): ^(x) can change in time by a phase 

d o M(x,y) = -> ^{x)=i){x + -x-)e~ lf{x++x ~\ (24) 

This phase can be chosen arbitrarily, without change of physics. We will choose a phase 
growing linearly with time f(x + + x~) = £(x + + x~)/2. Combining with (23) we arrive 
at 

This is the Hartree equation for relativistic particles interacting through a Coulomb poten- 
tial in the light cone variables. This can be seen more easily if we write (25) in momentum 
space. It will be more useful though to obtain (25) from a variational principle first. In 
fact (25) is equivalent to minimizing 

,2 

dxdy h(x — y)ip*(x)ip(y) 

(26) 

dxdyG{x-y)\^{x)\ 2 \^{y)\ 2 

under the constraint 

r dx\i;(x)\ 2 = 1. (27) 
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dx itfj* (x)d-ifj(x) — i—^- j i 



Notice that the funcional E is just the energy P + + P-, P+ being the hamiltonian H writen 
in terms of ifi. In momentum space E is given by 

,2 



(28) 



The symmetry ip(x) — > ip*(—x) of (25) implies that ip(p) is real for the ground state. Also, 
in order to have a positive definite energy ip(p) should vanish for negative p. This way we 
have 



4/ 



dpdqdk 1 r . . . . . . . . . , 2 / \2i 



(29) 
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-o [>p{k)ip{q)ip{k + p)ip{q-p) - ip{k) ip(q) ]. 



A discrete version of the energy above was minimized numerically using the method of 
the gradient [10], this means, starting from some arbitrary initial configuration the code 
changes it following the negative of the gradient of the discrete energy. The constraint 
(27) was imposed in two steps: 1) projecting out the component of the gradient normal 
to the constraint and 2) normalizing the state after each "time" step. Some examples of 
wave functions for diverse values of the parameter m/g are shown in figure 1. They are 
concentrated around a mean value of p- and are smooth, which implies that in position 
space ip vanishes at infinity faster than any polynomial. Notice that the two terms in E 
have competing effects. The kinetic term is minimized with a 5 function around p = m/2. 
The potential energy though favors a more spread out wave function. In fact we can see 
that the larger the value g for fixed m, the broader the function is. In the weak coupling, 
non relativistic limit m/g — > oo, the peak of the wave function is roughly around p = m/2. 
For the strong coupling region though, the peak is at larger values of p, reaching a value 
~ g/4 at m = 0. At this point the mass of the baryon is different from zero, signaling 
that the chiral limit does not correspond to m = (remember that m is the bare mass of 
the quark). This is disturbing and is probably related to the fact that in (21) the limit 
ty(x) = gives M(x, y) ^ 0. The ground state of the baryon number one sector is not only 
N c quarks on top of the vacum, but contains distortions even far away from the baryon. 
Notice it is the e(x — y) term in (21) that produces the mass renormalization back from \x to 
m. For negative values of m 2 the energy is still positive as long as m 2 is larger than some 
critical value ~ —2g 2 /3. Beyond this point E is not bounded from below. This behavior 
of the baryon mass as function of m/g is shown in figure 2. For large values of m/g the 
mass of the baryon (divided by N c , the number of quarks in the baryon) approaches the 
quark mass, as it should be for a non relativistic system. It is a little larger because in two 
dimensions the binding energy is positive. We also solved (23) linearised around the wave 
function ip of the (ground state) baryon. Now ip(p) does not need to be real (because 
it describes excited states) and we have one equation for the real part, another for the 
imaginary part. So, if 
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(30) 



then we have 

J j^G(k-p)(2Mk + q-p)Mq)*R(k) 

+ a R (k + q-p)Mq)Mk)) ( 31 ) 
= £a R {p) 

and 

(p+^W(p)-y / ^G(fc-p)(2^(fc + «-p)^(«y/(fc) 

- <7j(fc + g - p)ipo(q)ip (k)) ( 32 ) 

= 5(7/ (p) . 

The energies £ correspond to the excited states of the baryon. Typical examples of the 
wave function of the low lying states are shown in figure 3 and 4. Notice that there are 
two states nearly degenerate with the ground state. The source of this degeneracy is the 
presence of two symmetries in (25). The first is phase invariance 

%j)(x) -> e ia ^{x) (33) 

that implies that iipo(p) will be a solution with the same ground state energy. The second 
is translation invariance 

ip(x) — > ip(x — xq) (34) 

that implies that ipip(p) is also degenerate with the ground state. Translation symmetry 
is broken by the lattice so this degeneracy is slightly lifted in our numerical results, as can 
be seen in figure 3. The numerical procedure to find the n-th excited state was to minimise 
the energy in the subspace orthogonal to all states with lower energy, so instead of ipif^o (p) 
we have a linear combination of iptfjo(p) and iifJo(p) orthogonal to ipo(p)- 
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Figure Captions 

Figure 1 : Examples of wave functions of the ground state. P is measured in arbitrary 
units and g = 400. 

Figure 2 : Energy as a function of m/g for positive and negative values of m 2 (at the 
left of 0). Energy is measured in arbitrary units and g = 30. Negatives values energy are 
meaningless because the process of minimisation never reached an end. 
Figure 3 : Ground state and two other states degenerate with it. 
Figure 4 : Three examples of excited states. 
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